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Abstract. In the first part of this paper, we provide a concise review of our method of 
solution of the A^ Q-systems in terms of the partition function of paths on a weighted 
graph. In the second part, we show that it is possible to modify the graphs and transfer 
matrices so as to provide an explicit connection to the theory of planar networks introduced 
in the context of totally positive matrices by Fomin and Zelevinsky. As an illustration of 
the further generality of our method, we apply it to give a simple solution for the rank 2 
affine cluster algebras studied by Caldero and Zelevinsky. 
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1 Introduction 

Discrete dynamical systems may take the form of recursion relations over a discrete time variable, 
describing the evolution of relevant physical quantities. Within this framework, of particular 
interest are the discrete integrable recursive systems, for which there exist sufficiently many 
conservation laws or integrals of motion, so that their solutions can be expressed in terms 
of some initial data. Interesting examples of such systems of non-linear integrable recursion 
relations arise from matrix models used to generate random surfaces, in the form of discrete 
Toda-type equations [1, 17, 22]. More recently, a combinatorial study of intrinsic geometry 
in random surfaces has also yielded a variety of integrable recursion relations, also related to 
discrete spatial branching processes [4]. 

We claim that other fundamental examples are provided by the so-called Q-systems for Lie 
groups, introduced by Kirillov and Reshetikhin [20] as combinatorial tools for addressing the 
question of completeness of the Bethe ansatz states in the diagonalization of the Heisenberg spin 
chain based on an arbitrary Lie algebra. We proved integrability for these systems in the case 
of Ar in [7]. In the case of other Dynkin diagrams, evidence suggests integrability still holds. 

The Q-system with special (singular) initial conditions was originally introduced [20] as the 
recursion relation satisfied by the characters of special finite-dimensional modules of the Yan- 
gian Y{q), the so-called Kirillov-Reshetikhin modules. Remarkably, in the case q = Ar, the same 
recursion relation also appears in other contexts, such as Toda flows in Poisson geometry [15], 
preprojective algebras [14] and canonical bases [3]. 

In [7], we used methods from statistical mechanics to study the solutions of the Q-system 
associated with the Lie algebra Ar, for fixed but arbitrary initial conditions. Our approach 
starts with the explicit construction of the conserved quantities of the system, which appear as 
coefficients in a linear recursion relation satisfied by Q-system solutions. These are finally used 
to reformulate the solutions in terms of partition functions for weighted paths on graphs, the 
weights being entirely expressed in terms of the initial data. 

Note that there is a choice of various sets of 2r variables which constitute an initial condition 
fixing the solutions of the Q-system recursion relation, a choice parametrized by Motzkin paths of 
length r. This set of initial variables determines the graphs and weights which solve the problem. 
This is the key point addressed in [7], and is related to the formulation as a cluster algebra. 
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Cluster algebras [9] are another form of discrete dynamical systems. They describe a specific 
type of evolution, called mutation, of a set of variables or cluster seed. Mutations are rational, 
subtraction-free expressions. This type of structure has proved to be very universal, and arises 
in many different mathematical contexts, such as total positivity [13, 11], quiver categories [19], 
Teichmiiller space geometry [8], Somos-type sequences [12], etc. 

Cluster algebras have the property that any cluster variable is expressible as a Laurent 
polynomial of the variables in any other cluster in the algebra. It is conjectured that these 
Laurent polynomials have nonnegative coefficients [9] (the positivity conjecture). This property 
has only been proved in a few context-specific cases so far, such as finite type acyclic case [10], 
affine type acyclic case [2], or clusters arising from surfaces [24]. The Q-system solutions for Ar 
are also known to form a subset of the cluster variables in the cluster algebra introduced in [18] 
(a result later generalized to all simple Lie algebras in [6] ) . In [7] , we interpreted the solutions 
of the Ar Q-system in terms of partition functions of paths on graphs with positive weights: 
this proved positivity for the corresponding subset of clusters. Moreover, we obtained explicit 
expressions, in the form of finite continued fractions, for these cluster variables. We review our 
results in the first part of this paper. 

Of particular interest to us is the connection of cluster algebras to total positivity. Fomin 
and Zelevinsky [13] expressed a parametrization of totally positive matrices in terms of electrical 
networks, and established total positivity criteria based on relations between matrix minors, 
organized into a cluster algebra structure. In this paper we show the explicit connection of their 
construction to the Q-system solutions. 

In this paper, we review the methods and results of [7] in a more compact and hopefully 
accessible form. We first apply this method to the case of rank 2 affine cluster algebras, 
studied by Caldero and Zelevinsky [5]. These are the cluster algebras which arise from the 
Cartan matrices of affine Kac-Moody algebras. We obtain a simple explicit solution for the 
cluster variables in terms of initial data. We then proceed to describe the general solution of 
the Ar Q-system using the same methods. In particular, we obtain an explicit formula for 
the fundamental cluster variables, which generalizes the earlier results of [5] to higher rank. 
Finally, we make the explicit connection between the path interpretation of the solutions of the 
Q-system and a subclass of the totally positive matrices of [13] and their associated electrical 
networks. 

More precisely, it turns out that the generating function for the family of cluster variables Ri^n 
(n € N) of the Ar Q-system is given by the resolvent of the transfer matrix Tm associated with 
a graph Fm for some seed variables associated with the Motzkin path m. We show that it is 
possible to locally modify the graph without changing the path generating function, so that we 
obtain a transfer matrix of smaller size r + 1, equal to the rank of the algebra. From there, there 
is a straightforward identification with the networks associated with totally positive matrices of 
special type, related to the Coxeter double Bruhat cells of [15]. 

The advantage of our approach is that we have explicit expressions for the cluster variables 
in terms of any mutated cluster seed parametrized by a Motzkin path m. 

The paper is organized as follows. In Section 2, we explain the basic tools from statistical 
mechanics which we use, the partition functions of hard particles on graphs and the equivalent 
partition functions of paths on weighted dual graphs. 

For illustration, in Section 3, we use this to give the explicit expression for the generating 
function of the cluster variables of rank 2 cluster algebras corresponding to affine Dynkin dia- 
grams, a problem extensively studied in [27, 5, 25]. 

In Section 4, we review our solution [7] of the Ar Q-system for the simplest choice of initial 
variables. This solution uses the partition functions of Section 2 as well as the theorem of [23, 16] 
relating the partition function of non-intersecting paths to determinants of partition functions 
of paths. We also introduce a new notion in this section, the "compactification" of the graph. 
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which gives a new transfer matrix which is associated with the same partition function. This is 
a key tool in making the connection with totally positive matrices. 

Section 5 generalizes the results of Section 4, and we give expressions for the Q-system 
solutions in terms of any set of cluster variables in a fundamental domain. We also introduce 
the notion of "strongly non-intersecting paths" and a generalization of [23, 16]. This gives 
generating functions for the cluster variables in terms of the other seeds in the cluster algebra, 
and also provides a proof of the positivity conjecture [9] . This section is a quick review of the 
results of [7]. 

In Section 6, we extend the graph compactification procedure of Section 4 to the other 
graphs, corresponding to mutated cluster variables, introduced in Section 5. This yields transfer 
matrices of size (r + 1) x (r + 1), the resolvent of which is an alternative expression for the 
generating function for cluster variables. Finally, in Section 7, we use this to give the explicit 
relation to totally positive matrices [13, 15]. 

2 Partition functions 

2.1 Hard particles on Gr 

2.1.1 Vertex- w^eighted graphs 

Let G be a finite graph with N vertices labeled 1,2, . . . , N , and single, non-oriented edges 
connecting some vertices. The adjacency matrix A'^ of the graph G is the N x N matrix with 
entries Afj = 1 if vertex i is connected by an edge to a vertex j, and otherwise. To each 
vertex i is associated a positive weight y,. 

2.1.2 Configurations 

A configuration C of m hard particles on G is a subset of {1, . . . , A^} containing m elements, 
such that Afj = for all i,j £ C. 

This is called a hard particle configuration, because if we view the elements of C to be the 
vertices occupied by particles on the graph, the condition Afj = for i,j £ C enforces the rule 
that two neighboring sites cannot be occupied at the same time. Each vertex can be occupied 
by at most one particle. 

We denote by the set of all hard particle configurations on G with m particles. 

2.1.3 Partition function 

The weight of a configuration C is the product of all the weights associated with the elements 
of C. That is, 

The partition function for m hard particles on G is the sum over all configurations 
of the corresponding weights: 
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Figure 1. The graph G^, with 2r + 1 vertices labeled i = 1, 2, . . . , 2r + 1. 
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Figure 2. The graph G^, dual to G^, with total of 2r + 2 vertices. 



2.1.4 The graph Gr 

A basic example is the graph G = Gr of Fig. 1. It has 2r + 1 vertices and 3r — 1 edges, with the 
(symmetric) adjacency matrix defined by 

^1,2 = ^2r,2r+l = ^2j,2t+l = ^2j,2i+2 = ^2j+l,2i+2 = 1, iG{l,...,r — 1}. 

When r = 0, Go is reduced to a single vertex labeled 1, while for r = 1, Gi is a chain of three 
vertices 1, 2, 3 with the two edges (1, 2) and (2, 3). 

Example 2.1. The non- vanishing partition functions Z^'' for hard particles on Gr, r = 0, 1 are 



^0 



1, Zf°=yi, Z^^ = l, Z^-i = yi + 2/2 + ys, ^2^=2/12/3- 



7G1 



(2.1) 



2.1.5 Recursion relations for the partition function 



The transfer matrices .Z^'' satisfy recursion relations in the index r. They are obtained by 
considering the possible occupancies of the vertices 2r + 1 and 2r: 



(2.2) 



For example, Zq = 1 is the partition function of the empty configuration and Z^]^^ = \\ ?/2i+i 
for the maximally occupied configuration. 



2.2 Transfer matrix on the dual graph 

Associated to the graph Gr-, there is a dual graph Gr-, as in Fig. 2. 

It is dual in the sense that Gr is the medial graph of Gr- each edge of Gr corresponds to 
a vertex of G^, and any two edges of Gr share a vertex iff the corresponding vertices of Gr are 
adjacent. 

We fix the labeling so that the correspondence is between edges of Gr and vertices of Gr is: 

1) edge (/c, fc') of Gr corresponds to the vertex 2/c — 1 of Gr, where /c = 2, 3, . . . , r; 

2) edge (/c, A; + 1) of Gr corresponds to the vertex 2k of Gr-, where A; = 1, 2, . . . , r; 

3) edge (0, 1) of Gr corresponds to the vertex 1 of G^; 

4) edge (r + 1, r + 2) of G^. corresponds to the vertex 2r + 1 of G^. 
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r-1 



r+l r+2 




2' 3' 4' 5' f-l r' 

Figure 3. The graph Gr with oriented edges drawn in. 



2.2.1 Transfer matrix on Cjf'p 

We choose an ordering of the 2r + 2 vertices of Gj. to be < 1 < 2 < 2' < 3 < 3' < • • • < r < 
r'<r + l<r + 2, and consider this to be an index set. We construct the (2r + 2) x (2r + 2) 
transfer matrix Tr{ty) with these indices. Its entries are except for: 



{Tr)k',k = 1; {Tr)k,k' = *2/2fc-l) k 
{Tr)k+l,k = 1, {Tr)k,k+1 = ty2k, 

{Tr)l,0 = {Tr)r+2,r+l = 1, (Tr) 



0,1 



= 2,3,...,r; 
/c = l,2,...,r; 

{Tr)r+l,r+2 = ty2r+l- 



In matrix form, 



/O tyi 



Tr{ty) 




1 











\0 



ty2 


1 
1 












tys tyi 





































ty2r-l ty2r 
10 

10 ty2r+l 

1 / 



(2.3) 



This matrix is a weighted adjacency matrix for the graph Gr as drawn in Fig. 3. Each edge 
of Gr corresponds to two oriented edges pointing in opposite directions, with the weights in 
the transfer matrix corresponding to those oriented edges. The element of the transfer matrix 
indexed by i, j is the weight corresponding to the edge j ^ i. 

Lemma 2.1. The generating function for hard particle partition functions on Gr with weights 
tyi per particle at vertex i is 



r+l 



Z'^'-ity) := ^ t^'Zg^ = det(/ - T,(-ty)) 



(2.4) 



m=0 

with Tr as in equation (2.3). 



Proof. Expanding the determinant Dr{t) = det{I —Tr{—ty)) along the last column, one obtains 
the recursion relation: 

Dr{t) = Dr^l{t) + ty2r+lDr-l{t) + ty2rDr-2{t), r > 2. 

Let Z?r,m denote the coefficient of in Dr{t). It satisfies the same recursion relation as 
equation (2.2) for Z^--, with DqCO = 1 + ^2/1 and Di(t) = 1 + t{yi + y2 + ys) + f^yWs, in 
agreement with the initial values of Z'^^{ty), r = 0, 1 in equation (2.1). Thus, Dr^m = for 
all r, m. ■ 
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Example 2.2. For the case r = 1, Gi is a chain of 3 vertices 1, 2, 3. The dual Gi is a chain of 
three edges connecting four vertices 0, 1, 2, 3, and the transfer matrix is: 



nty) = 












1 





ty2 








1 







Vo 





1 


0/ 



One checks that det(/ — T{—ty)) = 1 + t{yi + y2 + ys) + i^yiy^, in agreement with the partition 
functions Zf'^ of equation (2.1). 



2.3 Hard particles and paths 

Let G be a graph with oriented edges. For example, the graph Gr of the previous section can 
be made into an oriented graph, by taking each edge to be a doubly-oriented edge. Then each 
oriented edge from i to j receives a weight w{i, j), which is the corresponding entry of the transfer 
matrix T(ty) of equation (2.3). This may be interpreted as the transfer matrix for paths on Gr 
as follows. 



2.3.1 The partition function of paths 

A path of length n on an oriented graph G is a sequence of vertices oi G, P = (uq, fi, . . . , Vn)-, 
such that there exists an arrow from Vi to Uj+i for each i. Let J'^(,(n.) be the set of all distinct 
paths of length n, starting at vertex a and ending at vertex h on the graph G. 

A path on a graph with weighted edges has a total weight which is the product of the weights 
w{vi,Vi+i) associated with the edges traversed in the path. 

The partition function Z'^^{n) for ^^^^{n) is 

n-l 

Assuming G is finite, and labeling its vertices i = 1, 2, . . . , A^, let us introduce the N x N transfer 
matrix T with entries Tij = w{j,i), we have the following simple expression for Zf^^{n): 

^Un) = (r"),_^. 

The matrix T(ty) of (2.3) is then the transfer matrix for paths on Gr with weights 1 for edges 
pointing away from the vertex 0, and tyt for the i-th edge pointing to the origin. The partition 
function for weighted paths of arbitrary length on Gr from the vertex to itself is 

^^^{ty) = E (^(*y)")o,o = ((^ - ^(*y))"')o,o- (2-5) 

n>0 

Lemma 2.2. The partition function of paths of arbitrary length from to on the graph Gr is 
equal to 

Proof. By definition, (2.5) is the ratio of the (0, 0)-minor of the matrix I — T{ty) to its determi- 
nant. Using equation (2.4) and the explicit form (2.3), this immediately yields the relation. ■ 

This relation between hard-particle partition functions and path partition functions may be 
interpreted as a boson- fermion correspondence, and is a particular case of Viennot's theory of 
heaps of pieces [28, 21]. 
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Figure 4. A path on Ga with 16 steps. The weights y,; are associated to the descending steps (1, —1) and 
to the second half of the horizontal steps (1, 0) + (1, 0) of the path, i being the label of the corresponding 
edge of G3 (see on left). Here, the path receives the weight y\y2VzV^yhV^V'j ■ 



2.3.2 The path partition function as a continued fraction 

A direct way to compute Z^''{ty) in equation (2.5) is by using Gaussian elimination on I — T{ty) 
to bring it to lower-triangular form. The resulting pivot in the first row is 1/ ((I — T(ty))~^)g ^. 
If we do this systematically, by left-multiplication by upper-triangular elementary matrices, ttie 
result is 

Lemma 2.3. 

^^'•(ty) = ^ ^ \, (2.7) 



2.3.3 Non-Intersecting paths: the LGV formula 

We may represent paths of length 2n from vertex to itself on Gr as paths on the lattice Z?,q (see 
Fig. 4 for an illustration). Such paths start at (0,0) and end at (2n, 0), and have the following 
possible steps: 

1) to the northeast, {j, k) — )• {j + 1, k + 1), corresponding to the jth. step in the path going 
from vertex k to vertex k + 1; 

2) to the southeast, (j, k) — )• (j + l,k — 1), corresponding to the jth step in the path going 
from vertex k to vertex k — 1; 

3) to the east, (j, k) — (j -|- 2, k), corresponding to the steps /c — A;' — /c if /c G {2, 3, . . . , r}. 

Paths from the origin to the origin have an even number of steps, by parity. 

We will need to consider the partition function of families of a non-intersecting paths on Gr, 
Zfj^. Here, the fixed starting points are parametrized by s = (si, . . . , So) and the endpoints by 
e = (ei, . . . ,ea)- 

By non-intersecting paths, we mean paths that do not share any vertex. We have the cele- 
brated Lindstrom-Gessel-Viennot formula [23, 16] 

2.?:= det 2,^'-„.. (2.8) 

The determinant has the effect of subtracting the contributions from paths that do intersect. 
This formula can be proved by expanding the determinant as: 

a _ 

5^sgn(a)nzf;„,e.. (2.9) 
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One then considers the involution ip on families of paths, defined by interchanging the beginnings 
of the two first paths that share a vertex until the vertex, or by the identity if no two paths 
in the family intersect. It is clear that when ip does not act as the identity it relates two path 
configurations with opposite weights in the expansion of the determinant, as the two starting 
points are switched by (p, hence these cancel out of the expansion (2.9). We are thus left only 
with non-intersecting families, all corresponding to cr = Id, hence all with positive weights. 



3 Application to rank 2 cluster algebras of affine type 
3.1 Rank two cluster algebras 

In this section, we use the partition functions introduced in the previous section to the prob- 
lem of computing the cluster variables of rank two cluster algebras of affine type with trivial 
coefficients. This allows us to give explicit, manifestly positive formulas for the variables, prov- 
ing the positivity of the variables in these cases. 



3.1.1 The recursion relations 

The rank 2 cluster algebras of affine type [9, 27] with trivial coefficients may be reduced to the 
following recursion relations for n G Z: 

1 



Cn+l 



-(1 + x^) if n is odd, 
+ if n is even. 



"7' (3.1) 



where 6, c are two positive integers with be = 4, hence (6, c) = (2, 2), (4, 1) or (1,4). The aim is 
to find an expression for n E Z in terms of some initial data, e.g. (xq, xi). 

2 -b^ 



The connection to rank 2 affine Lie algebras is via the Cartan matrix ^ 

~ C Zi 

The cases (4,1) and (1,4) are almost equivalent: If x„(xo,xi) is a solution of the (6, c) 
equation, then xi_„(xi, xq) is the solution of the (c, b) equation. We may thus restrict ourselves 
to the (1, 4) case. 

However, the symmetry Xn ^ a^i-n changes the parity of n. Therefore we need to also 
consider the dependence of Xn in the "odd" initial data (xi,X2). 

It turns out that the recursion relations (3.1) are all integrable evolutions. This allows us 
to compute the generating function for x„, n > 0. The result is a manifestly positive (finite) 
continued fraction. In the light of the results of the previous section, this allows to reinterpret Xn 
as the partition function for weighted paths on certain graphs. This path formulation gives yet 
another direct combinatorial interpretation for the expression of Xn as a positive Laurent polyno- 
mial of the initial data, to be compared with the approach of [27, 5] using quiver representations 
and that of [25] using matchings of different kinds of graphs. 



3.2 The (2, 2) case: solution and path interpretation 

Consider the recursion relation 

Xn+lXn~l = xl + 1 (3.2) 

with xo = (xo,xi) = (x,y). This is the Ai case of the renormalized Q-system considered in [6]. 
Due to the symmetry n o 1 — n of the equation, the solution satisfies 

x„(xo,xi) = Xi_„(xi,Xo) 

so we may restrict our attention to computing x„, for n > 0. 
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3.2.1 Constants of the motion 

Equation (3.2) is integrable. To see this, we rewrite (3.2) as 



2^n— 1 



1. 



Then 



= fn- 



We conclude that there exists a constant c, independent of n, such that Xn-i + Xn+i = cxn for 
all n. Using equation (3.2), 



Xn-l + Xn+l Xn+1 , 1 , X„ Xi 1 Xq 

c = = \ 1 = 1 1 . 

Xfi Xfi XnXn+l Xn+l Xq XqX\ X\ 



(3.3) 



We interpret c as an integral of motion of the three-term relation (3.2): all solutions of the latter 
indeed satisfy the two-term recursion relation (3.3), for some "integration constant" c fixed by 
the initial data. 

G 

Note that c coincides with the partition function Z^^ = yi + y2 + ys for one hard particle on 
the graph Gi, with weights 



yi 



xo' 



2/2 



XqXi 



ys 



Xq 
xi' 



(3.4) 



Note also that the only other non- vanishing hard particle partition functions on Gi are Zq ^ = 1 
and = yiys = 1. 



3.2.2 Generating function for Xn 

oo 

Let X{t) = ^ t^Xn be the generating function for the variables x„ with n > 0. Using Xn+i — 

n=0 

cXn + Xn-l = 0, we have by direct calculation 

Xq - {cXo-Xi)t Xq 



with yi as in (3.4). This gives Xn as a manifestly positive Laurent polynomial of {xq,xi). In 
fact, expanding the r.h.s. of (3.5), we get 



p,q,e>0 

from which x„ is obtained by extracting the coefficient of i". This agrees with the result of [5]. 

It follows that the dependence of the variables x„ on any other pair {xk,Xk+i) is also as 
a positive Laurent polynomial. This is clear from the translational invariance of the system: 



Xn{xk,Xk+l) = Xn--k{xo,Xl)\ 



x-oi-^a;fe,a;ii-4-a::fe_|_i 
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3.2.3 Relation to the partition function of paths 

Upon comparing the continued fraction expression (3.5) and equation (2.7), we see that there is 
a path interpretation to the variables follows. 

The denominator of the fraction X{t) is equal to the partition function Z'^^^—ty) for hard 
particles on the graph Gi introduced in Section 2, with the weights —tyi per particle at vertex i. 
Therefore, 

^f,. _ Z'^^{0,-ty2,-tys) 

° ZOi{-ty,,-ty2,-tysy 

so XQ^X(t) is equal the partition function for paths on the graph Gi of Example 2.2, beginning 
and ending at the vertex 0. The weights are as follows: The weight is equal to 1 for any step 
away from vertex 0, and and is equal to tyi per step from vertex i to i — 1. 

3.3 The (1,4) case: solution and relation to paths 

We now consider the system: 

1 + X2n-1 



X2n 



X2n-2 

1 + xi 



X2n+i = ^ ' ^" , neZ. 

X2n-l 

We determine the dependence of the variables Xn on two different types of initial conditions: 
case : (xq, xi), 

case 1 : (xi, X2). (3-6) 

We can eliminate the odd variables and get an equation for the even variables. Let Un — ^2n* 
Using X2n-i = UnUn-i — 1, the even variables satisfy the recursion relation 



Un+l 



+ Un-l 



UnUn-1 - 1 ' 

or, equivalently, ■u„('u„+in„_i - u^) = Un+i + Un-i- 

3.3.1 Conserved quantities 

The variable Wn = Un+iUn-i — u^ satisfies 

Un+l + Un-l u\ + Un-l + {UnUn-1 " '^)Un-l + /o ^7^ 

Wn = = 7 7T = -. (3.7) 

Un Un[UnU-i - 1) UnUn-l - 1 

Moreover, 

Un+l + Un+l + Un-l Un - Un+lUnUn-1 + Un+l + Un-l „ 

Wn+1 -Wn = = 7 TT = 0. 

Un+lUn - 1 Un Un[UnUn+l - 1) 

We conclude that Wn is a conserved quantity, that is, it is independent of n. 

Using (3.7), there exists a constant c such that Un+i + Un-i = cun- We may compute this 
constant explicitly in terms of the initial conditions in cases and 1 of (3.6), using uq = xq, 
ui = X2 and xqX2 = 1 + xi: 

n • .(0) - ^1+^0 _ ^o + ^i _ xl + {l + xif 

UiUq — 1 Xl XqXi 

(1) xj + jl + xif 
case 1 : ' = — ^ . 

XgXl 
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3.3.2 Generating function 

The linear recursion relation Un+i — cun + Un~i = implies the following formulas for the 
generating functions 

t/(»)(t) :,y„„f. = "°-"'''°'"°-"'' = ?2 (3.8) 



n>o 1 - tyi - — Tiy 

where the parameters are expressed in terms of the initial data as: 

(0) _ 1 + Xi (0) _ xl+{l + Xif (0) _ xl + l + xi 

Hi — 2 ^ y2 ~ 4 ' Us ~ 2 ' 

Xq XqX'^ XqX\ 

(1) _ Xj + I + Xl (1) _ xj + jl+Xif (1) _ 1+Xi 

ill — 2 ) y2 ~ 4 ' i/3 ~ 2 • l^O.iU; 

Note that both sets satisfy the same relation Ui y\ = 1 + where i = 0, 1, as they are 
a related via the substitution X2 -H- xq, which maps {y^i\y^\y^^^) ^ {y^^\y2'^ ^Vi'^) ■ 

Upon expanding the right hand sides of (3.8) and (3.9) in power series of t, we get all the 
UnS as explicit positive Laurent polynomials of either initial data (xo,xi) and {xi,X2)- Using 

we obtain the expressions, valid for all n > 0: 



l+i(q+e.)-2n-i{r+s) m-q-l 



case : X2n{xo,xi)= x^ " ' ' " ' 'x-^ 

q,£,r,s,m,>0 

- q - i\ /q + e - 1\ /n + 2r + s - 2q 
I — r,r J \ £ — s, s J \ m 



(3.11) 



easel: X2n+2{xi,X2) = ^ 



l+2n-i{q+l+r+s) q+t+m-n 



g,£,r,s,m>0 

+ + , + (3.12) 

\q — r,r, s J \ 1 J \ m 



where we have used the multinomial coefficients 



n 

mi, 771-2, • • • 



TTT-i! • • • mk\{n — ^ rrij)! 
otherwise 



The expressions for the odd variables follow from X2n+i = X2nX2n+2—^- The positivity of these 
latter expressions is easily checked: Note that the product X2n2;2n+2 in both equations (3.11) 
and (3.12) has a constant term 1 as a Laurent polynomial of (xo,xi) and (xi,X2) respectively, 
hence X2nX2n+2 — 1 is a positive Laurent polynomial. 
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(1.4) 



-m 



Figure 5. The graph 6(1^4) and the four corresponding path steps. 



3.3.3 Path interpretation 

The continued fraction expressions (3.8) and (3.9) allow for a path interpretation of X2n as 
follows. Consider the graph G(i 4) on the left hand side of Fig. 5, with two vertices labelled 0, 1 
connected by an edge, and connected to themselves via a loop. We assign weights to the oriented 
edges as follows: 



w{0 ^ 0) = tai 



w{0^1) = t, w{l^l)=ta3, w{1^0) = ta2. 



We can also associate a path in Z?,q to a path on G(i 4) composed of the steps shown in Fig. 5. 

The corresponding path transfer matrix is: T = t ) . Using Gaussian elimination on 

/ — T as in the previous example, we find that the partition function of paths from the vertex 
to itself is: 



((^-^no,o 



1 



1 — tai 



t^a2 
l-tas 



Hence, X2n is (up to a factor of xq in case and X2 in case 1) the partition function for paths 
of 2n steps on G(i 4), from and to the vertex 0, with weights (01,02,03) = (yi, 2/2,2/3) in case (i) 
and (01,02,03) = (^1,^2,-23) in case (ii), where the weights are as in (3.10). 



3.3.4 The graph Gi,4 as a compactification of the graph Gi 

We will see later that, quite generally, it is possible to "compactify" the graphs Gr- The result 
for r = 1 is precisely the graph Gi^4. Thus, The case (6, c) = (2, 2) of the previous section may 
also be interpreted in terms of paths on the graph G(i 4), but with different weights: 

w{0^0) = tyi, u;(0->l) = l, w{l ^ I) = tys, w{l ^ 0) = ty2. 

where yi are as in equation (3.4). 

To see this, note that the continued fraction (3.5) may be rewritten as: 

Thus, for n > 0, the solution Xn+i of the (2, 2) case is up to a factor xi the partition function 
for paths on ^(14), from and to the origin vertex 0, and with a total of n steps of the form 
0— )'0, 1— >-lorl— >-0, with respective weights yi, ?/2, 2/3. 
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4 Application to the Ar Q-system 

In [18, 6], we showed that the recursion relations (Q-systems) satisfied by the characters of the 
KiriUov-Reshetikhin modules of the quantum affine algebras associated with any simple Lie 
algebra can be described in terms of mutations of a cluster algebra. Solutions to the cluster 
algebra recursion relations are more general, in that the initial conditions are not specialized, as 
they are in the original Q-system satisfied by characters of Kirillov-Reshetikhin modules [20] . 

A particularly simple example of this is the case when the Lie algebra is of type A. Characters 
of An KR-modules, which are just Schur functions corresponding to rectangular Young diagrams, 
are given by the cluster variables upon specialization of the boundary conditions. In [7], we gave 
the general solution for the cluster variables without specialization of initial conditions. For 
these variables, we proved the positivity conjecture of Fomin and Zelevinsky by mapping the 
problem to a partition function of weighted paths on a graph. Let us review the results obtained 



4.1 Definition 

The Ar Q-system is the following system of recursion relations for a sequence Ra,n, a € Ir = 
{1, 2, . . . , r} and n G Z: 



with boundary conditions 

Ro,n = Rr+i,n = 1 for all n G Z. 

We note that the case r = 1 of coincides with the case (2, 2) treated in Section 3.2. 

Remark 4.1. The original Q-system, which is the one satisfied by the characters of KR-modules, 
differs from the system (4.1) not only in that it is specialized to the initial conditions Rafi = ±1 
for all a & Ir, but also by a minus sign in the second term on the right. In this discussion, we 
choose to renormalize the variables for simplicity. It is also possible (see the appendix of [6]) 
to consider the recursion relation with nontrivial coefficients. In that case both (4.1) and the 
original Q-system result from a specialization of the coefficients. 

We wish to study the solutions of equation (4.1) for any given initial data. Our standard initial 
data for the Q-system are the variables xq = {Ra,o, Ra,i)aeir- We may then view equation (4.1) 
as a three-term recursion relation in n, which requires two successive values of n (for each a) as 
initial data. Our first goal is to express all {Ra,n}a,n in terms of the initial data xq. 

4.1.1 Symmetries of the system 

We may use equation (4.1) to get both n > 2 and n < — 1 values in terms of xq. These 
are related via the manifest symmetry of the Q-system under n -H- 1 — n, which implies 
Ra,i-n{{Ra,o, Ra,i)a&ir) = Ra,n{{Ra,i, Rafi)a&ir) ■ In addition, equation (4.1) is translation- 
ally invariant under n ^ n + k: Ra,n+k{{Ra,k+i, Ra,k) o)aeir) foi' 3-11 



4.2 Conserved quantities of Q-systems as hard particle partition functions 



in [7]. 




(4.1) 



n,k £ 7j. 



The Q-system turns out to be a discrete integrable system, in that it is possible to find 
a sufficient number of integrals of the motion, as in the r = 1 case. 
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First, equation (4.1) may be used to express Ra,n (a > 2) as polynomials in {Ri^n}'- 

Ra,n = det (Ri^n+i+j-a-l)- (4.2) 

l<i,j<a 

This is proved by use of the Desnanot-Jacobi identity for the minors of the (a + 1) x (a + 1) 
matrix M with entries Mij = Ri^n+i+j-a, i,j £ la+i- 

The boundary condition i2j.+i,n = 1 (for all n) together with (4.2) implies the equation of 
motion which determines Ri^n- 



det {Ri^n+i+j-r-2) = 1- 

i,j<r+l 



From equation (4.1), it is clear that Rr+2,n = for all n. Hence there exists a linear recursion 
relation of the form 



r+l 



{-ircr+l-m Rl,n+m = 0, (4.3) 



m=0 



with Co = Cr+i = 1. The fact that c^+i-m do not depend on n follows from the fact that each 
row in the matrix in the determinant for Rf-^2 n is just a shift in n of any other row. 

The constants Cp,p = l,2,...,r are the r integrals of motion of the Q-system. They 
be expressed explicitly in terms of the -Ri,nS as follows: 



can 



Cp — det Rin+i+j~2 

' l<i<r + l I I I J 

l<j<r+2, j>^r+2-p 

independently of n. These quantities are similar to those found in [26] for the so-called Coxeter- 
Toda integrable systems. 

By using simple determinant identities, we show in Theorem 3.5 of [7] that Cp satisfy recursion 
relations which allow us to identify them as the partition functions Z^^ of p hard particles on 
the graph Gr of Fig. 1, where the weights given by: 

y2a-i,k = -B 5 , 1 < a < r + 1, y2a,k = — 5 5 , 1 <a<r. (4.4) 

^a,k'tla-l,k+l i^a,ki^a,k+l 

This is true for any k: The functions Z^^'" are independent of the choice of k. Thus, unless 
otherwise stated, will stand for yafi below. 

4.3 Q-system solutions and paths 

oo 

The linear recursion relation (4.3) allows to compute the generating function R^^^(t^ — t^R\^^ 

n=0 

explicitly. Indeed, R^''\t)Z^- {-ty) is a polynomial of degree r, and it is easy to see that 
Rir){t)Z^^{-ty) = i?i,oZ^'-(0, -ty2, -ty2r+i). 

Using (2.6), we may interpret R^^'\t) / Ri^ as the generating function 2,^Q(ty) for paths on Gr 
with the weights (4.4), say for k = 0. In other words, Ri^n/Ri,o is the partition function for 
weighted paths of 2n steps on Gr starting and ending at the origin (or starting at (0, 0) and 
ending at (2n,0) in the two-dimensional representation). 

Comparing the determinant formula for Ra,n (4.2) and the LGV formula (2.8) for families of 
a non-intersecting paths, we have [7] 
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/Cv^^\ 

Figure 6. The six pairs of non-intersecting paths on G2 of 8 and 4 steps, starting respectively at (0, 0) 
and (2,0) and ending at (6,0) and (8,0). 




1 2 ty^ 3 r ty.^^ r+1 



Figure 7. The graph GJ,, with r+1 vertices. We have indicated the weights attached to each oriented 
edge. 

Lemma 4.1. The quantity Ra,n/ {Ri.a)°^ is the partition function Z^^ for a non-intersecting 
weighted paths on Gr with starting and ending points Si = (2i — 2,0), = (2n + 2a — 2i,0), 
i = 1,2, . . . ,a. 

Proof. It is clear that Z'^.^^^. is the partition function for paths from {2i — 2, 0) to (2n+2a — 2j, 0), 
which is equal to that from (0, 0) to (2(n + a + 1 — i — j), 0) by translational invariance. This 
allows to identify the two determinants (4.2) and (2.8), up to an overall factor of (i?i^o)") and 
the result follows. ■ 

As an illustration, we have represented in Fig. 6 the six pairs of non-intersecting paths 
contributing to R2,3, solution of the A2 Q-system. 
Thus, we have proved 

Theorem 4.1. The variables Ra,n which satisfy (4.1), when expressed in terms of the va- 
riables xq, are equal to (Ri^q)" times partition functions for paths on the graph Gr, involving 
only the weights ya of (4.4) . These weights are explicit Laurent monomials of the initial data 
Xq = {Rafli Ra,i)aGir- This givcs an explicit expression for Ra,n cls Laurent polynomials of the 
initial data, with non-negative integer coefficients. 

4.4 An alternative path formulation 

In the same spirit as Section 3.3.4, one can show (see below) that the solution Ri^n of the A,. 
Q-system may also be interpreted in terms of paths on a new ( "compactif led" ) graph G'^ of 
Fig. 7. This is a graph with r + 1 vertices, labelled 1, 2, . . . , r + 1, connected via oriented edges 
i— )-i + l, i + 1— i = l,2, ...,r and loops i — )• z, i = 1, 2, . . . , r + 1. 
To each edge is attached a weight as follows: 

wt(z — ;> i + 1) = 1, wt(i + 1 -)• i) = ty2i, wt{i i) = ty2j-i- (4.5) 



The corresponding transfer matrix encoding these weights is an (r + 1) x (r + l)-matrix of the 
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form 



(tyi ty2 

1 tys tyi 
1 ty^ tye 



T' 



\ 



1 ty2r-3 ty2r-2 







Vo 

Then we have 



1 ty2r-l ty2r 

1 ty2r+ll 



R 



l + t 



yi 



1,0 



l-tyi- t- 



y2 



tys-t 



l-tyg-i- 



l-ii/2r-l-t 



52? 



(4.6) 



(4.7) 



This is readily proved by Gaussian elimination. 

Comparing this expression to the results of the previous section, we conclude that 

Lemma 4.2. Forn > 0, Ri^n+i/ is the partition function for paths on the weighted graph G'^, 
with a total of n steps along the edges of type i ^ i or i + 1 — )• i in G'^, starting and ending at 
vertex 1. 

The weights are given in equation (4.5) with yi as in (4.4) with k = 0. 

In Section 7 we will relate this result to the total positivity conjecture of Fomin and Zelevinsky 
and networks. 

We may actually write an explicit expression for Ri^n by simply expanding the continued 
fraction (4.7) as: 



(t) = Ri,o [l + tyi YI Yliiy2i+iy''^' 

,P2.--->P2r+l>0 e=0 
P0=P2r+2=0 

^P2i + P2e+1 + P21+2 - 1' 



X {ty2e+2r'+' 



P2£ - '^,P2i 



+1 



Substituting the values (4.4) for the weights ya = ya,o, and extracting the coefficient of t""*""^, 
we get for all n > 0: 



Rl,n-\ 



R 



1,1 



V- yr {R,^o)P^*+^+P^*+^-P^i-P^^-^ " fp2e + P2i+i + P2e+2 - 1 

Z-^ 11 ( -,\P2i+l+P2i-P2i-l-P2i-2 11 



Pl.P2.-'--P2r+l>0 i- 
P0=P2r+2=0. ^Pi=n 



i (i?- -^)P2i+l+P2i-P2i-l-P2i- 



i=0 



P2e - 1,P2<?+1 



as explicit positive Laurent polynomials of the initial data. This gives a rank-r generalization 
of the formula given in [5] for r = 1. 
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5 Cluster algebra formulation: mutations and paths 
for the Ar Q-system 

In this section, we show that the solutions {Ra,n \ ct G Ir,n G of the Q-system are positive 
Laurent polynomials when expressed as functions of an arbitrary set initial conditions. This 
generalizes our result for the initial condition xq in the previous section. 

The recursion relation (4.1) has a solution once a certain set of initial conditions is specified, 
but this set need not necessarily be the set xq. We will explain below that the most general 
possible choice of initial conditions is specified by a Motzkin path of length r. 

The solutions of (4.1) can be viewed as cluster variables in the Ar Q-system cluster algebra 
defined in [18]. Hence, our proof provides a general confirmation of the conjecture of [9] for this 
particular cluster algebra: When the cluster variables are expressed as functions of the variables 
in any other cluster (i.e. an arbitrary set of initial conditions), they are Laurent polynomials 
with non-negative coefficients. 

The results of this sections were explained in detailed in [7], and this section should serve as 
a summary of the proofs contained therein. 

5.1 The Q-system as cluster algebra 

In [18], it was shown that the Ar Q-system solutions {Ra,n} may be viewed as a subset of 
the cluster variables of the Ar Q-system cluster algebra. This is a cluster algebra with trivial 
coefficients, which includes the the seed cluster variable {Rifl, R2,o, ■ ■ ■ , Rr,o, ^2,1, ■ ■ ■ , Rr,i), 

with an associated associated 2r x 2r exchange matrix has the block form: , C the 

Cart an matrix of Ar- 

In this language, each cluster is a vector with 2r variables, and the subset of clusters relevant 
to the Q-system are those which have entries made up entirely of solutions to the Q-system 
(we restrict to these in the following). These clusters are all related by sequences of cluster 
mutations which are one of the relations (4.1). 

Because of the form of equation (4.1), it is easy to see that the restricted set of clusters 
corresponding to the Q-system are characterized by a set of r integers (mi, . . . ,mr), subject to 
the condition that \ma — ma+i\ < 1- This defines what is known as a Motzkin path. 

Definition 5.1. The cluster x^ corresponding to the set of integers m = (mi, . . . ,mr) is the 
vector of 2r variables {Ra,maj Ra,ma+i}aeirj ordered so that all variables with an even second 
index appear first. 

For example, the initial cluster xq corresponds to the Motzkin path m = (0, 0, . . . , 0), and 

xo = {Rlfl, • • • , Rrfl, Ri,i, ■ ■ ■ , Rr,l)- 

For any m, Xm is obtained from the fundamental initial seed xq by TfiutatioTis of the cluster 
algebra. This is just saying that one gets Xm by repeated selected applications of the recursion 
relation (4.1) to xq. Each mutation changes only one of the cluster variables. That is, for some a 
and n, 



Recall that we only consider here particular mutations that only involve solutions of the Q- 
system. Here we have used the "time" variable n to define forward (resp. backward) mutations 
according to whether the mutation increases (resp. decreases) the index n in the mutated cluster 
variable. 




(forward mutation at a), 
(backward mutation at a). 
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a 







(a) 



(b) 



(c) 



Figure 8. The Motzkin path m = (2, 1, 2, 2, 2, 1, 0, 0, 1) for r = 9 (a) is decomposed into p = 6 descending 
segments (12)(3)(4)(567)(8)(9) (circled, red edges). The corresponding graph pieces are indicated 
in (b) . They are to be glued "horizontally" for flat transitions (green edges) and "vertically" for ascending 
ones (blue edges). The resulting graph is represented (c) with its vertex (black) and edge (red) labels. 

Alternatively, the mutation changes the Motzkin path which characterizes Xm, by changing 
m I— )• m lb e^, with the plus (minus) sign for a forward (backward) mutation. Here, Ca is the 
vector which is zero except for the entry a, which is equal to 1. We see that accordingly the 
Motzkin path locally moves forward (backward). 

The Laurent property of cluster algebras ensures that every cluster variable is a Laurent 
polynomial of the cluster variables of any other cluster in the algebra. 

The positivity property, proved only in particular cases so far, is that these Laurent polyno- 
mials have non-negative integer coefficients. The property was proved for the particular clusters 
considered in the present case in [7]. It may be stated as follows: 

Theorem 5.1 ([7]). Each Ra,n, when expressed as a function of the seed Xm for any Motzkin 
path m, is a Laurent polynomial of \^Ra m^, 

Ra,mc+i I o £ Ir}, with non-ncgativc integer coeffi- 
cients. 

We outline the proof below. 

For clarity let us introduce the following notation. Let F be some cluster variable. Then F 
can be expressed as a function of Xm for any m. The functional form is then denoted by 
F = Fni(xin). Since F can also be expressed as a function of any other cluster, we can write 
-Pm(xm) = -^*m'(^m') for any two Motzkin paths m, m'. In particular, in the notation of the 



5.2 Target graphs and weights 

Due to the reflection and translation symmetries of the Q-system, we can restrict our attention 
to seeds associated with Motzkin paths in a fundamental domain = {m | m.ma{ma.) = 0}. 
There are 3^~^ elements in Mr- 

To each Motzkin path m E M,., we associate a pair (Lm, {ye(m)})) consisting of a rooted 
graph Tin with oriented edges, and edge weights ye(m) along the edges e. 
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5.2.1 Construction of the graph Fm 

The graph Fm is constructed via the following sequence of steps (see Fig. 8 for an illustration) : 

1. Decompose the Motzkin path m into maximal "descending segments" rtij of length ki 
{i = l,...,p). These are segments of the form rrij = {ma^,ma^+i, . . . ,ma^_^i-i) with 
Oi+i = Oi + ki, where ruai+j = ruai — j- Here, ai = 1 and a^+i — 1 = r. 

2. The separation between two consecutive descending segments of the Motzkin path, rrij and 
rrij+i is either "flat" i.e. rua^^-^ = ^0^+1-1 or "ascending" i.e. rna^_^_-^ = mQ,-_^j_i + 1. 

3. To each descending segment nij, associate a graph T^., which is the graph G^. with 
additional, down-pointing edges a ^ b for all a, b such that /i;j-|-l>a>5-l-l>2. There 
are a total of ki{ki — l)/2 extra oriented edges. 

4. We glue the graphs Fm; and Fmj+i into a graph FniJ|Fni,_|_i defined as follows (see Fig. 8 
for an illustration): 

(a) If the separation between nij and rtij+i is flat, we identify vertex of Fmi+i with 
vertex ki + 2 of T^., and vertex 1 of F„ii+i with vertex /c, + 1 of Fm. , while the 
connecting edges are identified. 

(b) If the separation is ascending, we reverse the role of vertices and 1 in the procedure 
above. 



The result of this procedure is the graph Fm = Fini||Fm2|| • • • ||rmp- Its root is the vertex 

of Fmi- 

We label the vertices of the graph Fm by the integers i with i S {0, . . . ,r + 2 + 714.(111)} 
(where n+(m) is the number of nia such that rria+i = rua + 1) and labels i' for any univalent 
vertex attached to vertex i via a horizontal edge. We do this by labeling the vertices of Fj^^ 
from bottom to top, by shifting the labels of the subgraphs Fm^ so that no label is skipped nor 
repeated. 

The edges e pointing towards the root of Fm are of two types: 

{i) the "skeleton edges" belonging to some G^. in the above construction; 
(n) the extra, down-pointing edges added in the gluing procedure. 

5.2.2 The weights on the graph Fm 

We label the 2r -|- 1 skeleton edges of type (i) by a = 1,2, . . . ,2r + 1 from bottom to top (see the 
example in Fig. 8), and the weights are denoted by ?/a(m). Weights assigned to edges pointing 
away from the root are all set to 1. 

Alternatively, we may label the "down pointing" skeleton edges by the pairs of vertices 
i + l^i or i'^i which they connect. The extra edges of type (ii) are also labeled by the pairs 
a — )• 6 of vertices which they connect. All edge weights may be labeled by the label of the edge. 

The weights of the edges of type (ii) can be expressed in terms of the skeleton weights: 




n yi+i,i(m) 



b<i<a 



b<i<a 



SO that they obey the following intertwining condition 



ya,6(m)ya',6'(m) = Vafi' {m)ya' ,b{m) , 



a > a > b > b' . 



(5.1) 
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For example, the extra weights of the example of Fig. 8 read respectively: 2/3^1 = ^2^4/2/8) 
y9,7 = yi4yi2/yi3, 2/8,6 = yuyio/yn, and yg^Q = ygjysfi/yu = yuyuyio/iyisyn)- 

Finally, for a given Motzkin path m G M^, we define the skeleton weights ?/ci,(m), a = 
1,2, . . . ,2r + 1 to be: 

y2a-i(m)= ^"'"^^ , a = l,2,...,r + l, (5.2) 

I N _ /^a+l.rfia + l / -1 c- Aa+i^mc+i r- 
y2almj — I i - dma^ma_+\+\ + -\ (>mc,m„+i+l 

X ( 1 - '^mc-i,m^+i + <5m^_i,mc+i ) , a = l,2,...,r, (5.3) 



where 



^ -^Q,ri+1 Ra,n 



Note that with these definitions the expressions (5.2), (5.3) involve only variables of the seed x^n. 

To each Motzkin path m G M,., we may finally associate a transfer matrix Tm = T^{ty{ra)), 
with entries {Tm)b,a = weight of the oriented edge o — )• 6 on Fm- Then the series in t 

Zmit;^m) ■■= ((/-r^)-i)oo = j;rzj'j(n) (5.4) 

n 

is the generating function for weighted paths on Fm, with the coefficient of t'^ being the partition 
function of walks from vertex to itself on Fm which have n down-pointing steps. When 
m = mo = 0, this coincides with (2.5). 

5.3 Mutations, paths and continued fraction rearrangements 

Our purpose is to write an explicit expression for the functional dependence of the variables Ra,n 
on the seed variable Xm, that is, find (i?a,n)m(xm) for each a, n and a Motzkin path m. 
To do this, we will describe how the generating function R^'^^(t^ — -^rn^(^;^m 

) is related to 

the generating function i?[^/(t; x^')) where m and m' are related by a mutation. Then we start 

(r) / \ (r) 

from the known function -Rmo(^j ^mo)) ai^d apply mutations to obtain all other functions Rm 
with m in the fundamental domain Mr. 

One can cover the entire fundamental domain starting from iuq = by using only forward 
mutations m 1— )• m' = m + ea of either type (i) (. . . , a, a, o + 1, . . . ) 1— )• (. . . , a, a + 1, a + 1, . . . ) 
and (a) (. . . , a, o, a, . . . ) I— )• (. . . , o, a + 1, a, . . . ), with the obvious truncations when a = 1 or r. 
(See Remark 8.1 in [7].) 

Suppose m and m' are related by such a mutation. We compute the two generating functions 
of the type (5.4), Zm{t;Xm) and ZmK*! ^m')- 

In fact, the two matrices, Tm and T^' differ only locally, so that in computing the two 
generating functions by row reduction, we find that the calculation differs only in a finite 
number of steps. Note that generating functions take the form of finite continued fractions with 
manifestly positive series expansions of t. 

We note two simple rearrangement lemmas which can be used to relate finite continued 
fractions: 

{R2) a+- = - JJ-, wheve a =a + b, b= — — -, c= — — . 

1 — c 1 — 1-^ a + b a + b 

One checks this by explicit calculation. 
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Let a > 1. Then, using {R2) we can show that 
weights yi = yi{m) and y'^ = yi{m') are related via: 



2'm'(*;^m') if and only if the 



(ii) TUa-l = rua = rUa+l 



y2a-l+y2a, /3 = 2a - 1, 

y2ay2a+l/{y2a-l + y2a), (3 = 2a, 

y2a-iy2a+l/{y2a~l+y2a), /3 = 2a + 1, 

y^, otherwise, 

y2a-l+y2a, /? = 2a - 1, 

y2ay2a+l/{y2a-l + y2a), P = 2a, 

y2a-iy2a+l/iy2a-l+y2a), /3 = 2a + 1, 

2/20+22/20- l/(y2Q-l + y2a), /3 = 2a + 2, 

y^, otherwise. 



One checks directly that the expressions (5.2), (5.3) indeed satisfy the above relations. 

The boundary case, where a = 1, is treated analogously, but first requires a "rerooting" of 
the graph to its vertex 1, which is implemented by the application of (i?i): Indeed, we simply 
write Zm{t;x^) = {{I - Tm)'^)^^ = 1 + tyi(m)Z,'^(t; x^) with Z'^{t;x^) = {{I - T^)'^)^^^. 
We then rearrange 2.'^(t;xm) using {R2) again, and find that Z'j^{t;Xm) = 2,m'(^; ^m') if and 
only if the weights are related via the above equations. 

The net result is the following. Given a compound mutation /^m which maps the fundamental 
Motzkin path mo to m = {mo,)a£i^, then there are exactly nii "rerootings" as described above. 
This corresponds to rewriting the generating function 



-Km 



mi — 1 

(i;Xni) = ^ fRi^i + t"^^Ri^rniZm{t;^r. 
1=0 



with 2.m(i;Xni) as in equation (5.4). 

This leads to the following main result: 

Theorem 5.2. For each n > 0, the function (i?i.n+mi)m(xm) = Ri,mi Z'o'o (^)- Thus it is 
proportional to the generating function for weighted paths on the graph Fm with positive weights, 
so it is a manifestly positive Laurent polynomial of the initial data x„i- 



We have represented in Fig. 9 the graphs Fm for the Motzkin paths m of the fundamental 
domain for r = 3. 



5.4 Q-system solutions as strongly non-intersecting paths 

To treat the case of Ra,n with a 7^ 1, given a Motzkin path m G M.,., we need a path interpre- 
tation for the determinant formula for Ra n+mi '■ 

i?o,n+mi _det(i2l,n+mi+.+,-o-l)i<,,<, /^^^ _ ^^^^^ 



{Rl,nii)°' Rl,n+mi ^ ' ^ ^<i,j<a 

Here, we have used the result of the previous section to rewrite the formula in terms of the 
partition function for paths on Fm, from and to the root, and with n + i+ j — a — 1 down steps. 
As in the standard LGV formula, we interpret this determinant as a certain partition function 
for paths on F„i starting from the root at times 0, 2, 4, . . . , 2a — 2 and ending at the origin at 
times 2n, 2n + 2, . . . , 2n + 2a - 2. 
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Figure 9. The Motzkin paths m of the fundamental domain M3 and the associated graphs Fm, with 
their vertex and edge labels. We have also indicated the mutations by arrows, the label being a when 
the mutation fi^ acts on variables Ra,m with an even index m and a + r for an odd index m. 

5.4.1 Paths on rGprGSGntcd cis paths on a square lattice 

We draw paths, with allowed steps dictated by the graph Fm, on a square-lattice in two dimen- 
sions. Paths start and end at y-coordinate 0. Moreover, if a path has n "down" steps (steps 
towards the vertex 0), then its starting and ending point are separated by 2n horizontal steps. 
That is, a path is from (x, 0) to (x + 2n, 0) where x the starting time and n is the number of 
down-steps. 

Since the horizontal distance between the starting and ending points is fixed by the number 
of "down" steps, a single step of the form 

a — )■ 6 = a — /i 

should have a horizontal displacement 2 — h (instead of 1 as in the usual case). That is, on the 
square lattice it is a segment of the form 

(x, a) — )• (x + 2 — /i, a — h). 

Some examples are illustrated in Fig. 10. 

Thus, we identify Z-qq" (n + — a — 1) = 2,^5^^ with paths on the two-dimensional lattice 
starting at the point Sj = {2i — 2, 0) and ending at the point ej = (2n + 2a — 2j, 0), with the 
types of steps allowed given by the edges of Fm in the way explained in the previous paragraph. 
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target 



up steps 



down steps 




path 



Figure 10. The two-dimensional representation of a typical path on the graph Fm, m the strictly 
descending Motzkin path (2,1,0) of the case ^3. Descents of /i = 2 are vertical (time displacement by 
2 — ft, = 0), while descents of /i = 3 go back one step in time (time displacement by 2 — ft = —1). With 
these choices, the total time distance between start and end is twice the number of descents (16 = 2x8 
here) . 



5.4.2 Strongly non-intersecting paths 

We now look at families of a paths, corresponding to the determinant in equation (5.5). Such 
paths may have crossing on the lattice. As in the case of LGV formula, the determinant cancels 
out contributions from paths which share a vertex. However, other situations may occur: Two 
paths may cross without sharing a vertex in our picture. 

One can generalize the proof of the LGV formula to take such crossings into account. Using 
the expansion (2.9), and introducing an involution on families of paths. This involution 
interchanges the beginnings of the first two paths which share a vertex or which cross each 
other, by transforming the crossing segments [P, Q\ and [i?, S\ into non-crossing ones [i?, Q\ 
and [P, 5] . This effectively interchanges the two paths up to the points P and R respectively, 
whichever comes first. 

As in the usual case, the involution 99 acts as the identity if no two paths cross, share a vertex, 
or can be made to cross via such an exchange. 

Taking into account the weights of the paths, the intertwining condition (5.1) implies that 
the flip preserves the absolute value of the weight, but changes its sign, due to the transposition 
of starting points. So the determinant (5.5) cancels not only the paths that share a vertex or 
that cross, but also those that come "too close" to one-another, namely that can be made to 
cross via a flip. 

We call the families of paths which are invariant under the involution ip strongly non- 
intersecting paths. 

To summarize, we have the following theorem: 

Theorem 5.3. For any Motzkin path m G M^, the variable Ran+mi {with n + mi > a — 1) 
is equal to (-Ri,mi)" times the partition function of a strongly non-intersecting paths with steps 
and weights determined by T^. The starting points are Si = (2i — 2,0) and the end points are 
Cj = (2n + 2a — 2j, 0), with i,j = l,...,a, and the weights are functions of the cluster x„i. 
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(a) (b) (c) 

Figure 11. The identification of horizontal (a) or vertical (b) pairs of consecutive vertices on Fm and the 
result on F'^. The new edge on F^ with weight —1 allows to subtract the contribution from paths that 
do not exist on Fm- We have also represented the situation of a longer chain (c), where the identification 
now requires a network of up-pointing edges with alternating weights ±1 for the suitable subtractions. 

In particular, (i?a,n+mi)m(xm) is a Laurent polynomial with non-negative coefficients of the 
cluster Xin. 

6 A new path formulation for the Ar Q-system 

In Section 5, we constructed a set of transfer matrices T^, associated with paths on the 
graphs Fm, which allowed us to interpret Ri^n as generating functions of weighted paths on 
a graph Fm, and hence prove their positivity as a function of the seed variables Xm. 

In Section 4.4, we also showed that for the special case m = mo = 0, there is an alternative 
graph F'^p = G'^, and that one can interpret Ri^n as a generating function for paths this 
alternative "compactified" graph. The graph F^^ has r + 1 vertices, hence the associated 
transfer matrix is of size r + 1 x r + 1. 

We now ask the question, is there a corresponding compactified set of graphs, F^, which 
give a path formulation of Ri^n with weights which given by functions of the seed Xm for all 
Motzkin paths in M^.? 

It turns out that it is always possible to find a weighted graph with r + 1 vertices which 
answers this question positively for each m. This corresponds to a set of transfer matrices of 
a size equal to the rank of the algebra Ar. Therefore this transfer matrix allows us to make 
a direct connection between our transfer matrix approach and the totally positive matrices 
of [13]. 

6.1 Compactified graphs 

Consider the collection of graphs Fm. If we are interested in the generating function of weighted 
paths on them from the vertex to itself, then we can make various changes in them locally 
( "compactify" them) without affecting the generating function itself. 

We obtain such compactified graphs from Fm by identifying pairs of neighboring vertices, 
and, when necessary, adding oriented edges to cancel unwanted terms. There are two possible 
ways to make such identifications. Before presenting the general case, let us illustrate the two 
situations in the following subsection. 
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6.2 Examples of compactification 

In fact, the first type of compactification, applied to Gr, leads to the alternative graph G'^ 
obtained in Section 4.4. Recall that the generating function for paths from 1 to 1 on the 
graph Gr of Fig. 3 is related to the generating function for paths from to via the rerooting 
procedure: 

Example 6.1. The generating function \ is equal to the generating function of paths from 1 
to 1 on the graph obtained from G^ via the following procedure ("compactification"): 

1. Identifying vertex i' with vertex i, whenever both exist, and attaching a loop with weight 
tyii i to the resulting vertex (see Fig. 11 (a)). 

2. Identifying vertex r + 2 with vertex r + 1, and attaching a loop with weight tyr.+2,r+i to 
the resulting vertex. 

3. Identifying vertex 1 with vertex 0, renaming the resulting vertex 1, and attaching a loop 
with weight ty\ = tyi,o to this vertex. 

This is clear, since a path from 1 to 1 through the vertex i' must have a segment i — )• i' — )• z with 
weight ti/i'^i, and if it goes through the vertex r + 2 it must have a segment r + 1 — 7'r + 2— )-r + l 
with weight tyr+2,r+i- Similarly, the loop at 1 accounts for segments of the form 1 — )• — )• 1. 
Note that the resulting graph is G'^ of Fig. 7. Thus, we have the identity of Section 4.4 

^1,1 — ^1,1- 

Another example of identification of vertices is the case when the two vertices are adjacent 
vertices, + 1), on the spine of F (see Fig. 11 (b)). 

Example 6.2. Consider the graph associated to an ascending Motzkin path segment of 
length k (see the example for k = 3 m the lower right hand corner of Fig. 9). This is a vertical 
chain of 2A; + 2 vertices, numbered from to 2/i; + 1 from bottom to top, connected by edges 
oriented in both directions. The edges i + 1 — )• z have weights yj+i and the edges i ^ i + 1 have 
weights 1. 

Suppose we identify the vertices a, a + 1 in this graph for some a. A path from to with 
a step a — )• a + 1 is always paired with a step a + 1 — )• a, for a net contribution to the weight 
of the path of i/a+i- We associate a loop with weight i/a+i at the newly formed vertex after the 
identification of vertices a and a + 1. 

However, there are "forbidden" paths on the resulting graph, paths which are not inherited 
from paths on Hr- These are paths which go from a+2 — t- a — 1 without traversing the loop. This 
would correspond to going from a + 2 to a — 1 on Hr without passing through the edge a + 1 — t- a, 
which is impossible. We cancel the contribution of these paths by adding an ascending edge 
a — 1 — )• a + 2 with weight —1. The effect is precisely to subtract the weights of the forbidden 
set of paths. 

More generally, a succession of identifications of the type in the example above results in the 
following (see Fig. 11 (c)): 

Lemma 6.1. The generating function of paths from to on the graph Hj. with vertices 
0, . . . , 2A; + 1 is equal to the generating function for paths from to on the following com- 
pactified graph H'j^ : 
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1. Identify the vertices 2i + l and 2i + 2; rename the resulting vertex i + 1 (i = 0, . . . ,k — 1). 

2. Attach a loop at the vertex i + 1 with weight ifj+i = y2i+2- Other edge weights remain 
unchanged. 

3. Add ascending edges j — )• j + 2 + a (0 < a < — 1 — j, < j < A; — 1) with weight (—1)''"'"-^ 
to the resulting graph. 

An illustration of the resulting graph H'j^ is given in (c) of Fig. 11. 

Proof. Consider the set of paths 7 of the form PiP^ P2P~ P^, where Pi are fixed paths, and P"*" 
is a path from /iq to hi consisting of only up steps and loop steps, and P^ is a path from hi 
to /iQ consisting only of down steps and loop steps. We furthermore restrict ourselves to paths 
with the weight(P^)x weight(P~) fixed to be yw" = w^^ ■ ■ ■w^'^ for some n. Here, y is the 
product of the weights of the down steps in P~ , and w" is the total weight coming from the 
loops in the path. Let / be the weight of the remaining fixed portions of the path. 
Without loss of generality, we can take /iq = and hi = k + 1. 

For each such path we can decompose Ui = nf + , where nf is the number of times the 
loop with weight Wi is traversed in P"*" , and in P~ . Paths which arise from paths on Hj. must 
have > 1 (for all i) by definition. We claim that on H'f^, paths with n~ = are cancelled by 
paths which pass through the new ascending edges. 

The key observation is that a path which has an up step going through the ascending oriented 
edge i — 1— )-i + a + l (a>0)on iJ^ has nf = n^-^ = . . . = m+a = 0. 

Then the total contribution of the paths in T to the partition function is in fact 

fy E w" = /y E ^" 

-{-n . —rii nT^+n.— ^l,■ 

^ i * I I <■ 

n~>0 n+>0 



fy E w"-/yE E 



nT=0 
3 



+ /yE E -"-/y E E + 

jl <j2 n+ =n+ =0 ii <h <h n+ =n+ =n+ =0 

Jl 32 n 32 J3 

= /yE(-ir E E 

a=0 h<--<ja n+ =■■■=«+ =0 

31 3a 

So the alternating sum has the effect of subtracting the terms with any n~ = 0. 

A path on the graph i/^ with a spine vertices skipped (by traversing ascending edges of length 
> 1) comes with a total sign (—1)^ where x = ^(length of the ascending segments — 1) = a. 
That is also the sign of the term with nf = • • • = nt = in the summation. 

Jl Ja 

Finally, we note that any path can be decomposed into pairs of ascending and descending 
segments as above, and the proof can be applied iteratively to any path. ■ 



6.3 Definition of compactified graphs 

On a graph Fm, we call a skeleton edge horizontal if it connects (a) vertices i and i' for some i, 
(b) vertices and 1, or (c) the top vertex and the one below it. We call an edge ueriimZ otherwise. 

Definition 6.1. The compactified graph F^ with r + 1 vertices is obtained from the graph Fm 
via the following compactification procedure: 
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(a) (b) (c) 

Figure 12. (a) The graph Fm of Fig. 8, and (b) the compacted graph F^. We have circled on Fm the 
pairs of vertices to be identified in the compactification procedure (in red for horizontal pairs, in blue 
for vertical pairs). We also represent in (c) the Motzkin path m. 

1. Introduce an order on the vertices of Fm, so that i < i + I and i < i' < i + 1. Number 
them from 1 to 2r + 2 accordingly. 

2. Identify vertices 2j — 1, 2j (j = 1, . . . , r + 1), and rename the resulting vertex j. Double 
edges connecting (2j — 1, 2j) are replaced by a loop at j with weight which is the product 
of the weights on the two edges. All other edges and their weights are unchanged. 

3. All maximal subgraphs of the form Hk of Fm, consisting of vertical edges only, are replaced 
by compactified weighted graphs of the form ff^, as in Lemma 6.1 (with the obvious shift 
in labels). 

Example 6.3. For illustration, the identification of edges in the case of the graph of Fig. 8 (c) 
are: ~ 1, 2 ~ 2', 3 ~ 4, 5 ~ 5', 6 ~ 6', 7 ~ 7', 8j- 8', 9 ~ 9', 10 ~ 11 and 12 ~ 13. See Fig. 12. 
There are two maximal subgraphs of the form Hi are the vertices 2, 3, 4, 5 and 9, 10, 11, 12. 
Each pair now corresponds to a vertex i in F'j^, i = 1, 2, . . . , 10, which receives a loop i — )• i from 
the identification. 

Fig. 13 shows the set of compactified graphs corresponding to Fig. 9 for the case A3. 
The resulting weighted graph F'^ has r + 1 vertices labelled l,2,...,r + l, and hence is 
associated with a transfer matrix of size r + 1 x r + 1. 

6.4 An alternative construction 

An alternative description of the compactified graphs F^ is the following. 

We start from the graph F'^^ = G', of Fig. 7. The loop at vertex i has weight ty2i-i and the 
edge z + 1 — )• i has weight ty2i, where i/j = yj(mo) are as in equations (5.2), (5.3). 

Decompose m into maximal segments of the form: 

1) descending segments, Da,i = ((a, m), (a + 1, m — 1), . . . , (q + i — 1, m — i + 1)); 

2) ascending segments, Aa^i = {{a, m), (a + 1, m + 1), . . . , (a + i — 1, m + i — 1)); 

3) flat segments ((a, m), {a + 1, m) , . . . , (a + k — 1, m)). 

Here, i > 2 and k > 1. 
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Figure 13. The Motzkin paths m of the fundamental domain M3 and the associated compact graphs F^, 
with their vertex and edge labels. Note the up-step weights: yi^s — ?/2,4 = ^1 and yi_4 ~ 1. Mutations 
are indicated by arrows. 

Definition 6.2. The graph T'^ is the graph obtained from F^^j = G'^ via the following steps: 

1. For each descending sequence Da^i we add descending edges Q+p— )-a + g(0<g + l< 
p < i) to F'j^jj, with weights tya+p,a+q{'ioa.), where 



ya+p,o+g(m) 



a+p—1 

n y2i(m) 

j=Q+9 

a+p—1 

n y2i-i(m) 

j=a+q+l 



(6.1) 



2. For each ascending sequence Aa^i, we add ascending edges a+q — )■ a+p (0 < q + 1 < p < i) 
with weights ya+q,a+pim) = (-1)^"^"^ 

Lemma 6.2. The weighted graph F'^ is identical to the weighted graph F'^. 

This is just the result of the definition of Fm using the decomposition of m, as in Fig. 8. 
Maximal subgraphs the form correspond to the maximal ascending segments of the Motzkin 
path. All other segments correspond to subgraphs with horizontal edges. 



6.5 Equality of generating functions 

To summarize, the compactified graph F'^^^ is such that 
Theorem 6.1. 

(i-T^)rj = (i-Orj- 

In other words: the partition function for weighted paths from vertex 1 to vertex 1 in Fm is 
identical to that for weighted paths from vertex 1 to vertex 1 in the compact graph F'^. 
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7 Totally positive matrices and compact if ied transfer matrices 

We now establish the connection between the transfer matrices for paths on T'^ and the 
totally positive matrices of [13] corresponding to double Bruhat cells for pairs of Coxeter ele- 
ments. 

We may express the compact transfer matrices of the previous section in terms of the 
elementary matrices fi, e^, di for GL^+i, defined as follows. Let Eij denote the standard 
elementary matrix of size (r + l) x (r + l), with entries {Eij)k.i = 5k,i^i,j- 

Definition 7.1. The elementary matrices {ej,/j,(ij} are defined by 
/i = / + Ai£'i+i,i, ei = I + ViEi^i+i, i G {1, . . . ,r}, 

di = I + - l)E,^i, i G {l,...,r + l}. (7.1) 
for some real parameters Aj, /ij, Ui. 

In [13], Fomin and Zelevinsky introduced a parametrization of totally positive matrices as 

r+l 

products of the form 11 /« H '^i H J suitable sets of indices, and A,, /ij, Vi some 

ie/ i=i jej 

positive parameters. This expression allowed to rephrase total positivity in terms of networks. 
Here we interpret our compact transfer matrices in terms of some of these products. 

Recall that each Motzkin path can be decomposed into descending, ascending and flat pieces, 
as in Section 6.4. We introduce the increasing sequence of integers (oi, . . . , a2fc)i such that the jth 
ascending piece of m, of m starts at rria- = a2j-i and ends at ma^+ij-i = a2j- Similarly 

for the sequence of increasing integers (bi, . . . , 62^), which mark the starting and ending points 
of the descending sequences Da^,i- 

For i < j, let w[i,j] denote the permutation which reverses the order of all consecutive 
elements between z and J in a given sequence. That is, a;[i, j] = (i, j)(i + 1, j — l)(i + 2, j — 2) • • • . 
For example, uj[i,j]: - 1, . . . , i) (i, i + 1, . . . , j), 




(Jr 



-1,...,1). (7.2) 

Example 7.1. For the Motzkin path m = (2,1,2, 2, 2, 1, 0, 0, 1) of Fig. 8, we have the ascending 
segments [2,3] and [8,9], while the descending segments are [1,2] and [5,7]. The rearranged 
sequences read a^n = (8, 9, 7, 6, 5, 4, 2, 3, 1) and = (9, 8, 5, 6, 7, 4, 3, 1, 2). 

Note that the sequences cJm and Tm consist of increasing and decreasing subsequences of 
consecutive integers, and that these subsequences and their order are unique. 

One can define the decomposition of the transfer matrix into a strictly lower-triangular 
part and an upper triangular part B^, so that 

= + ^m- 

Lemma 7.1. The matrices o-nd can be expressed as 

Nm = I-{fnfi,---fi^r\ (7.3) 
Bm = t{did2 ■ ■ ■ dr+i) (ejiCja • • • Cj,.), (7.4) 

where the parameters in equation (7.1) are 

Xi = 1, fii = 1/21-1, J^i = (7.5) 

2/2 j-l 
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Figure 14. Pictorial representation of the action of (a) f~^^ (b) an increasing product fa^fa+i ' ' ' fb-i^ 
and (c) a decreasing product f^\f^\ - ■ ■ fa^ for some a < b. These correspond to adding (a) an 
ascending edge j — J' j + 1 with weight —1, (b) ascending edges a + i ^ a + i + 1, i = 0,1, . . . ,b — a — 1 
with weights —1, and (c) a web of ascending edges a + i^a + k,0<i<k — l<b — a — 1 with weights 



Proof. We give a pictorial proof. It is possible to describe multiplication by an elementary 

matrix as the addition of an arrow to a graph. In our context, A'^m encodes the ascending arrows 

in the graph F^, and Bm the descending arrows. 

First, consider the product f^^ ■ ■ ■ f^^ in Am of equation (7.3), where {fi)~^ = I — 

The sequence (v, V-i, • • • 1 ^1) = ( Hi ^[0-21-1: CL2i\) o (1, . . . , r), which is simply written in 

reverse order, consists of alternating increasing and decreasing sequences of consecutive integers. 

The products of matrices corresponding to increasing subsequences are 

a2j+i— a2j — 1 
1=1 

and products corresponding to decreasing sequences are 

a2j—a2j-i i—l 

Pj — fa2jfa2j-l" ' fa2j-i ~ ^ ~^ ^ -E'a2j_l+j,a2j_i+fe- 

i=l k=0 

We start with the graph corresponding to the identity matrix, which is the transfer matrix 
of the graph consisting of r + 1 disconnected vertices labelled l,2,...,r + l, each with a loop of 
weight 1. Multiplying on the left by /~ creates an ascending edge j — ?• i + 1 with weight — 1. 
More generally, left multiplication by creates a succession of ascending edges 02^+1 — a2j+2, 
a2j + 2 — )• a2j +3, . . . , 02^+1 — 1 — )• a2j+i- Left multiplication by P~ creates a web of ascending 
edges a2j-i + fc — )• a2j-i +i, 0<k<i — 1< a2j — a2j-i — 1, with alternating weights (—1)*+'^. 
We illustrate the resulting actions on the graphs in Fig. 14. 

Recall that the segments [a2i~i,a2i] correspond to the ascending segments of m, themselves 
associated to the vertical chain-like pieces of Fm (see Fig. 8). Recall that in the identification 
procedure leading to F^ (Definition 6.2 and Lemma 6.2), we showed that such chains must 
receive a web of ascending edges with alternating weights ±1 (see Fig. 11 (c)), while all the 
vertices are connected via ascending edges i — t- i + 1 with weight 1. 

Finally, comparing this with the graph associated to (/ii/jj ' " fik)~^ described above, 
we find that the contribution of ascending edges to (or equivalently, F'j^) is identical to 

^ ~ {fiifi2 ' ' ' fik) 
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The proof of equation (7.4) is similar, but now concerns the descending edges and the loops 
of the graph T'^. 

The product tdid2 ■ ■ -c^r+i is the transfer matrix of a chain of r + 1 disconnected vertices 
i = l,2,...,r + l, each with a loop with weight ty2i-i. Multiplication on the right by creates 
a descending edge i + 1 — >• i, with weight ty2i-i = ty2i- 

y2z — 1 

Again we divide the sequence Tm into increasing subsequences of consecutive integers, (62^-1, 
62J-1 + 1, • ■ • , ^2j) and decreasing subsequences (^2^+1 — 1, • • • , ^2j + !)• Therefore the product 
ej-^ • • • ejr consists of "ascending" factors 

b2j — b2j~i i—1 

Qj = e-b2j-l(^b2j ■ ■ ■ eb2j = 1+ ^ ^ 2/b2^._j+j,f,2j_i+fc^62j_l+fc,62j~l+i 

i=l k=0 

and "descending" factors 

b2j+i—b2j — l 

„_ r , y2b2j+2i 

— efe2j + l-l6&2j + l-2 • • • eb2j + l — J + > iib2 .+j,fe2j+j+l) 

~[ y2b2j+2i-l 

where yb+ifi+j are the weights of equation (6.1). 

Recall that the segments [62j-i)^2i] correspond to the descending segments of m, which 
correspond to networks of descending edges on Tm with weights (6.1) (see Fig. 8). In our 
construction 6.2 of F^, these descending edges have remained unchanged, while each vertex i 
received a loop with weight ty2i-i, « = l,2,...,r + l. This is nothing but the graph associated to 
t{did2 • • • dr+i){ej-^ej^ ■ ■ ■ ej^), which therefore encodes the contribution of loops and descending 
edges to T^, and equation (7.4) follows. ■ 

We can now make a direct connection with totally positive matrices encoding the networks 
associated to the Coxeter double Bruhat cells considered in [15]. Each Motzkin path m G 
corresponds to such an element, the (r + 1) x (r + 1) matrix P^: 

Definition 7.2. Given a Motzkin path m G Mr, define 
-Pm = ifiifi2 ■■■fir) {did2 ■ ■ ■ dr+l) (ejiCja • • • e^J, 

where (1112 ■ ■ - h) = fm and {jij2 ■ ■ ■ jr) = Tm are the two sequences of (7.2). The parameters 
Xj, Hj, Uj are as in equation (7.5). 

As a consequence of equations (7.3), (7.4), we have 

Theorem 7.1. 



(I - = {i- tp^r' (/. J.2 • • • /v), 

which allows to rewrite the generating function of Theorem 6.1 

'((/-tPm)-^)l^l+ ((/-ti^m)-^)l_2 if 

02+1 

{(I-tPm)-'),^, if ai = l. 



a=ai 



This yields an interpretation of the solution Ri^n+mi+i to the Ar Q-system with initial 
data Xm in terms of the network associated to the totally positive matrix for all n > 0. 
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Figure 15. The network corresponding to the totally positive matrix Pm for the Motzkin path m of 
Example 7.1, depicted in Fig. 12 (c). We have indicated in medallions the three network representations 
(in red, black, blue) for the elementary matrices fi, di and e^. The network for corresponds to the 
same picture, but with the / part (red descending elements) to the right of the e part (ascending blue 
elements). 



Example 7.2. For the fundamental Motzkin path m = mo with ma = for all a, we have 
o"mo ="^1110 = - 1, • • • , 1), and therefore P^o = (/r/r-i • • • /i)(dic^2 • • • 4+i)(erer-i • • • ei). 
Note that the matrix F = frfr-i ■ ■ ■ fi has entries Fij = 1 if i > j, and otherwise. One can 
check directly that I — tPmo = ~ ^mo)' with T^^^ given by equation (4.6). 

An equivalent formulation uses the explicit decomposition Pm = FDE, where F = fi^ - ■ ■ fi^, 
D = di - ■ ■ dr+i and E = ej^ - ■ ■ ej^ : 

CO oo 

(I - tP^y^F = f'iFDEYF = fF{DEFY = F{I - tP^y\ 

n=0 n=0 

where P^ = DEF = p-^P^F, and the fact that F is a lower uni-triangular matrix, which 
implies: ((I - tT^Y^)^^^ = ((/ - tP^)-^F)^^^ = ((/ - tP^)-^)^^^. 

The matrix P!^ is another way to write a totally positive matrix, and the network graph 
corresponding to it has a slightly modified form from that of P^n- Both of these correspond to 
electrical networks [13]. For illustration, we represent in Fig. 15 the network corresponding to the 
matrix P^ for the Motzkin path m of Example 7.1, represented in Fig. 12 (c). The medallions 
show the three elementary circuit representations for the three types of elementary matrices 
di , 6^ , each receiving the associated weight. The Lindstrom lemma [23] of network theory states 
that the minor |P|ri|'.!'.|rt of the matrix P of the network, corresponding to a specific choice rows 
ri, . . . , rfc and columns ci, . . . , Cfc, is the partition function of k non-intersecting (vertex-disjoint) 
paths starting at points ri, . . . , and ending at points ci, . . . , c^, and with steps taken only on 
horizontal lines or along f,dore type elements. Here we have only considered circuits with 
entry and exit point 1 , after possibly several iterations of the same network (each receiving the 
weight t), and whose generating function is precisely the resolvent ((/ — tP)~^)^ ^. 
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Figure 16. The concatenation of n copies of the network coded by P^. The quantity (P^i)" i the 
partition function of electrical wires starting and ending at the two indicated arrows. We also represent 
below a cylinder formulation a la [15]: the wire must wind n times around the cylinder before exiting. 



8 Conclusion 



In this paper, we have made the contact between our earlier study of the solutions of the A^- 
Q-system, expressed in terms of initial data coded via Motzkin paths, and the totally positive 
matrices for Coxeter double Bruhat cells. We showed in particular how the relevant pairs of 
Coxeter elements were encoded in the Motzkin paths as well. 

One would expect the total positivity of the transfer matrices Pm or P!^ to be directly related 
to the proof of the positivity conjecture in the case of the Q-system. Our proof presented 
in [7] relies on the path formulation of and on the formulation of Ra,n as the partition 
function of families of strongly non-intersecting paths. The total positivity of the compactified 
formulation should provide an alternative proof, using networks rather than paths. 

The precise connection between paths on graphs and networks, as illustrated in Section 7 
above is subtle. Indeed, the identity between resolvents implies that the partition function for 
weighted paths from 1 to 1 on Tm with n descents, is identical to the generating function 

for circuits on a network made of n identical concatenated networks, each corresponding to the 
totally positive matrix P^, from connector 1 to connector 1 (see the top of Fig. 16). In [15], 
this concatenation is realized by putting the network on a cylinder and allowing for the circuit 
to wind n times around it before exiting (see the bottom of Fig. 16). Note that we could also 
work with P^ instead, as it is related to P[^ via cyclic symmetry. 

More generally, it should be possible to relate our non-intersecting path families to networks 
with multiple entries and exits, as in the setting of the Lindstrom lemma. 

Another question concerns the cluster algebra attached to the Ar Q-system. As stressed in [7], 
we have only considered a subset of the clusters which arise in the full Q-system cluster algebra, 
namely those which consist of solutions of the Q-system. There are other cluster mutations, 
however, which are not recursion relations of the form (4.1). One may ask about the other cluster 
variables in the algebra. The positivity conjecture should hold for them as well. Preliminary 
investigations show that the corresponding mutations can still be understood in terms of (finite) 
continued fraction rearrangements, hence we expect them to also have a network counterpart. 
These clearly can no longer correspond to Coxeter double Bruhat cells, as those are exhausted 
by the solutions of the Q-system. 

Finally, the connection to total positivity should be generalizable to the case of other simple 
Lie algebras as well. Indeed, on the one hand the Q-systems based on other Lie algebras also 
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have cluster algebra formulations [6], while on the other hand the notion of total positivity 
has been extended to arbitrary Lie groups [11]. We have evidence that hard particle and path 
interpretations exist for all Q-systems, and it would be interesting to investigate their relation 
to the corresponding generalized networks. The integrability of these systems is presumably 
related to that of the Coxeter-Toda systems of [26]. This will be the subject of forthcoming 
work. 
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